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The thermodynamic cycle of a gas parcel in the thermo-acoustic engine is referred to as a thermo-acous¬ 
tic micro-cycle, which consists of two isobaric branches by two straight line branches. The influence of 
quantum degeneracy on the work output of the cycle is investigated based on the correction equation 
of state of an ideal Bose gas. The relationship between the dimensionless work output W* and the effi¬ 
ciency rj* is obtained under the condition of weak gas degeneracy. Some significant results are discussed. 

© 2008 Elsevier Ltd. All rights reserved. 


1. Introduction 

A thermo-acoustic engine (prime mover and refrigerator) [1-4] 
is of the advantages of high reliability, low noise, simple construc¬ 
tion, non-parts of motion, non-pollution, ability to self-start, etc. It 
can utilize a wide variety of energy resources: solar, geothermal, 
industrial waste heat and marsh gas. It has very important signif¬ 
icance for environmental protection and moderating the tense 
petroleum needs in the world. With this great potential, the ther¬ 
mo-acoustic engine has captivated many scientists and engineers 
in the thermal science and power and cryogenic engineering. 

Since few decades, a lot of research works, initialed by Rott [5], 
dealing with thermo-acoustic, render the growing interest in 
studying interaction phenomena between thermal and acoustical 
waves. There are predominantly more standing-wave thermo¬ 
acoustic engines involved in wide research. The idea of a travel¬ 
ing-wave thermo-acoustic engine was put forward more than 20 
years ago [6]. In general case, the acoustical filed will be in be¬ 
tween a pure traveling-wave and a pure standing-wave. It is note¬ 
worthy that traveling-wave and standing-wave modes are quite 
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different. There are different thermodynamic processes and differ¬ 
ent expressions for acoustical power production, thermal efficiency 
for the traveling-wave acoustic system and the standing-wave 
acoustic system [7,8]. 

We assume that the working gases consist of many gas parcels. 
A single parcel oscillates close to its equilibrium point by means of 
acoustic wave. The universal oscillation characteristics for each 
parcel of fluid are alike. The dissimilarity between the traveling- 
wave mode and the standing-wave mode is depicted by unlike 
phase displacement between the presses oscillation and the veloc¬ 
ity oscillation. 

In the 20th century, several researchers had believed that the 
thermodynamic cycles of the gas parcels in a thermo-acoustic en¬ 
gine consist of reversible adiabatic steps and irreversible constant- 
pressure steps, which is identical to the Brayton cycle [1,2]. But the 
analysis shows that the parcel has always contacted with the plate 
as it oscillates alongside the plate. The adiabatic steps are not tally 
with the fact. In this paper, a new cycle model will be adopted to 
study the thermo-acoustic cycle. 

According to the theory of classical thermodynamics, the 
performance parameters of the cycle can be derived by using the 
classical ideal gas equation of state or some equations based on 
it. However, when the gas density is high enough or temperature 
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is low enough, the ideal gas will deviate from classical gas behav¬ 
iour and quantum degeneracy of the gas will become important 
[9-17]. Under the conditions, the gas called quantum gas obeys 
Bose-Einstein statistics instead of classical statistics. One can call 
the cycle with quantum Bose gas as quantum cycle or Bose cycle. 
Thus, when the working fluid is a quantum gas, the cycle will have 
some new performance characteristics, which are different from 
those with conventional working fluid. 

The objective of this paper is to establish the cyclic model of a 
typical gas parcel at the middle of the regenerator in a thermo¬ 
acoustic engine. Unlike a reversible Carnot cycle, the performance 
of a thermo-acoustic micro-cycle is dependent on the properties of 
the working fluid. Under the high initial pressure the working fluid 
may become ideal quantum gases, harmonic oscillator systems, 
spin systems, and so on. This paper analyzes the effect of quantum 
degeneracy on the performance of the thermo-acoustic cycle 
which consists of two constant-pressure branches by two straight 
line branches. Some special cases are discussed in detail. 
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2. Thermo-acoustic cycle of a gas parcel at the middle of the 
regenerator 

The working fluid of a thermo-acoustic engine consists of a lots 
of gas parcels [1,18]. The thermodynamic cycle of a gas parcel in 
the thermo-acoustic engine is called as thermo-acoustic micro-cy¬ 
cle. Consider a parcel of gas at the middle of the regenerator that 
moves back and forth along the plates at the acoustic frequency. 
As it moves, the parcel of gas will experience changes in tempera¬ 
ture. The presses oscillation in a real thermo-acoustic engine is 
sinusoidal. According to the general model of thermo-acoustic cy¬ 
cle [1,2], we can instead of the sinusoidal wave with the correction 
square- wave for simplicity as shown in Fig. 1. The constant pres¬ 
sure may be regarded as an average value between state 1 and 
state 2 of the sinusoidal oscillation. Figs. 2 and 3 show the ther¬ 
mo-acoustic micro-cycle of a gas parcel. When the parcel moves 
from the high temperature end with temperature T x + VT of 
the plate to the state 2 with temperature T 2 = T 0 + T m alongside 
the plate (see Fig. 2a), the heat passes into the parcel due to the 
temperature of the plate being higher than that of the parcel. The 
volume of the parcel increases and its pressure is simplified as a 
constant p 0 + p m . When the parcel moves from the state 2 to the 
low temperature end with temperature T x - £ m VT of the plate 
(see Fig. 2a), the heat puts into the plate due to the temperature 



Fig. 2. Thermodynamic processes of a gas parcel in the thermo-acoustic engine. 

P 



Fig. 3. Pressure-volume diagram of the thermo-acoustic micro-cycle. 


of the plate being lower than that of the parcel. The volume of 
the parcel remains increase and its pressure decreases. When the 
parcel moves from the low temperature end to state 4 with tem¬ 
perature T 4 = T 0 - T m (see Fig. 2b), the heat is put into the plate 
due to the temperature of the plate being lower than that of the 
parcel. The volume of the parcel decreases and its pressure is sim¬ 
plified as a constant p 0 - p m . When the parcel moves from the state 
4 to the high temperature end of the plate (see Fig. 2b), the heat 
passes into the parcel due to the temperature of the plate being 
higher than that of the parcel. The volume of the parcel continues 
decrease and its pressure increases. In Fig. 2, Sk is the thermal pen¬ 
etration depth, T 0 , Po, T m , p m and are the mean temperature, the 
initial pressure, the temperature amplitude, the pressure ampli¬ 
tude and the displacement amplitude of the parcel, T x and VT are 
the local temperature and the temperature gradient along the 
plate. 

For a sinusoidal wave, the presses can be written as [7,8] 

P(t) =P 0 -hPd sin (cot + 9 ) (1) 

where 9 is phase displacement between the presses oscillation and 
the velocity oscillation. For the correction square-wave oscillation, 
p m may be assumed as 

Pm = Pml cos 9 + p m2 sin 9 (2) 

where p ml and p m2 are constants. It is noted that 6 is a constant for 
given thermo-acoustic system, thereby p m is a constant. 
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The pressure-volume diagram of the thermo-acoustic cycle of a 
gas parcel is shown in Fig. 3. These processes 1-2 and 3-4 are iso- 
baric, 2'-3 and 4-T are isothermal, 2-Q-3 and 4-Q-l are adia¬ 
batic, 2-3 and 4-1 are straight line. The cycle l'-2'-3-4-l' is 
identical to the Ericsson cycle, the cycle 1-2-Q-3-4-Q-1 is identi¬ 
cal to the Brayton cycle, and the cycle 1-2-3-4-1 is the cycle mod¬ 
el built in this paper. 


3. The thermodynamic properties of the gas parcel 


The gas parcel consists of many Bose particles. These expres¬ 
sions for number density and pressure of ideal Bose gases are given 
as follows [19]: 


n N gf3/ 2 {z) 

V 2 3 

( 3 ) 
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( 4 ) 


with T 2 = T 0 + T m and T 4 = T 0 - T m . Using Eqs. (4) and (13), one can 
obtain 
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By using Eq. (10) for an isobaric process, the heat quantity Q u 
supplied by the plate in process 1-2 and the heat quantity Q 34 re¬ 
leased to the plate in process 3-4 are, respectively, 
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From Eq. (9) the change of internal energy in process 2-3 is given by 
A u 23 = jH C v dT = | Nk[T 3 F(T 3 ,p 3 ) - T 2 F(T 2 ,p 2 )] (17) 


with 


Using Eq. (11) yields the work in the process 2-3 as following: 
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where g is the number of possible spin orientations of a particle, k is 
Boltzmann’s constant, T is the temperature of gas parcel, 
X = {2nt £ kT)0 .5 is the mean thermal wavelength of the particles, h is 
Planck’s constant, //(z) is called the Bose function, z = exp {p/kT) is 
the fugacity of the gas, r(l) is the Gamma function, and p is the 
chemical potential. Using Eqs. (3) and (4) yields 

p = nkTF(z) ( 6 ) 


where F(z) can be called as the correction factor and it is given by 
[19] 
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with 


0(T - T e ) 


1 (T > T e ) 
0 (T < T e ) 



where T e is the Bose-Einstein condensation temperature. 

Based on quantum statistics and Eq. ( 6 ), the heat capacity at 
constant volume and the heat capacity at constant pressure can 
be expressed as 


c ’Aw)ri Nk w mT - V) ' 
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The equations described straight line processes 2-3 and 4-1 can 
be expressed as 
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with Pi = p 2 = Po + Prn and p 3 = p 4 = p 0 - p m . 

On the other hand, for an ideal Bose gas the ratio of the entropy 
to the total number of particles N is only a function of the ratio of 
the chemical potential /a to temperature. Because the states 2 and 3 
are at the isentropic line 2-Q-3 and the states 4 and 1 are at the 
isentropic line 4-Q-l, one can obtain 


F(T 2 ,p 2 ) = F(T 3 ,p 3 ), F(T,,p,) = F(F 4 ,p 4 ) (13) 


Combining Eqs. (13), (17) and (18) gives the heat quantity Q 23 
released to the plate in the process 2-3 as following: 


0.23 — I At/23 + W23 
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Combining Eqs. (9) and (12) with (13) and (17) gives the heat quan¬ 
tity Q 4 i supplied by the plate in the process 4-1 as following: 


Q 41 — Au 4 i + W 4 1 
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4. The work output and the efficiency of the cycle 


Combining Eqs. (13)-(16), (19), (20) yields the net work of the 
cycle W= Q .12 + 0.41 - IQ 23 + Q 34 I ^s follows: 
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The heat quantity Qi = Q .12 + Q 41 absorbed from the plate in a 
cycle is 
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From Eqs. (21) and (22), the efficiency of the cycle p = W/Ch can be 
obtained. 

Under the condition of weak gas degeneracy and first approxi¬ 
mation, the correction factor for Bose gas can be simplified as [20] 


F(T,p) = 1- 


Dp 

fW2 


(23) 


with D = /(4V2k 25 ), where h is Planck’s constant. 

Substituting Eq. (23) into Eqs. (21) and (22) yields 
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w * = 
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with x = bn. and y = ^. In Eqs. (24) and (25), one has neglected the 
third-order wave quantity. 


5. Discussion 


(1) It is noteworthy that the efficiency rf is independent on x in 
Eq. (25). Combining Eq. (24) with (25) gives the relationship 
between the dimensionless work output W* and the 
efficiency rf 


W* 




1.5 Dp f 
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The dimensionless work output W versus the efficiency r) 
characteristic with p 0 = 800,000 P a and T 0 = 20 K are shown in 
Fig. 4, where we have assumed the thermo-acoustic engine in 
space shuttle is at the condition of outer space. One can see 
clearly from the curves in Fig. 4 that the dimensionless work 
output W has a maximum value W* m . W* m increases with the 
increasing of x. Taking the derivatives of W with respect to p* 
and setting it equal zero, one can find that when p* = 
the dimensionless work output W* approaches maximum value 



0.625 



1.5Dp c 

t2.5 




It can be clearly seen from Eq. (27) that the maximum dimen¬ 
sionless work output W* ]m increases with the increase of x 2 . 

(2) The effects of the initial pressure p 0 on the dimensionless work 
output VV with T 0 = 20 1< and x = 0.05 are shown in Fig. 5. It is 
seen from the curves in Fig. 5 that the dimensionless work out¬ 
put W* has a maximum value W pm , too. Obviously, for different 
parameter p m , l/V versus p 0 characteristics be different, but 
these values at maximum dimensionless work output W pm 
are alike.In Eq. (24), taking the derivatives of Vf with respect 
to the initial pressure p 0 and setting it equal zero, one can find 
that when p 0 satisfies the following equation: 

Po = (1.25^-0.75^) ( 28 ) 

the dimensionless work output W* approaches maximum 
value 


w m = 0.225 fe) + 0.75^ + 0.625x 2 (29) 

Vo / 1 o 

The corresponding efficiency of the cycle can be expressed 
as 



0.25x-+!+ 

_ ‘o 
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It shows that the optimal dimensionless work output W pm is 
monotonically increasing function of the temperature ampli¬ 
tude T m and the pressure amplitude p m . 

(3) Under the limit 4^ -> 0, the correction factor goes to zero, Eq. 
(26) gives 


W* 


5 ( * \ ( x _ 2 F ^ 

VI - 2 ri*)\ \-2vr) 




Fig. 4. Dimensionless work output W versus the efficiency r\ with p 0 = 800,000 P a 
and T 0 = 400 K. Curve 1: x = 0.05, Curve 2: x = 0.06, Curve 3: x = 0.07. 



Fig. 5. Dimensionless work output W versus the initial pressure p 0 with T 0 = 400 K 
and x = 0.05. Curve 1: p m = 30,000 P a , Curve 2: p m = 35,000 P a , Curve 3: 
p m = 40,000 P a . 


In this case, the cycle becomes a classical thermo-acoustic cycle. 
Combining Eq. (26) with (31) gives that the dimensionless work 
output Vf is larger than that of the classical ideal gas. 


6. Conclusion 

The thermodynamic model of a thermo-acoustic micro-cycle for 
a typical parcel at the middle of the regenerator has been estab¬ 
lished in this paper and the optimal dimensionless work output 
of the cycle using an ideal Bose gas as work fluid has been ana¬ 
lyzed. The following important conclusions are obtained. 

(1) A thermo-acoustic micro-cycle can be simplified as a model 
which is made of two isobaric branches by two straight line 
branches. The straight line branches tally with the fact better 






















F. Wu et al.f Applied Energy 86 (2009) 1119-1123 


1123 


than the adiabatic branches. The heat exchanging quantities, 
the work and the change of internal energy of the straight 
line processes are calculated under the condition of weak 
gas degeneracy. 

(2) The initial pressure p 0 should be selected according to Eq. 
(28) to obtain the maximum dimensionless work output of 
the thermo-acoustic micro-cycle. 

(3) There exist maximum values of the dimensionless work out¬ 
put W* related to the efficiency p* and the initial pressure p 0 . 
The analysis shows that the maximum dimensionless work 
output W pm is monotonically increasing function of the tem¬ 
perature amplitude T m and the pressure amplitude p m . 

(4) Under the condition of weak gas degeneracy, the efficiency 
of the cycle p* is the same as that of the classical ideal gas, 
while the dimensionless work output Vf is larger than that 
of the classical ideal gas. 
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